Breathers In Periodic Granular Chains With Multiple Band Gaps 
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We consider the localized nonlinear breathing modes that emerge in heterogeneous granular con- 
figurations of two materials with a periodicity of three and four beads. We examine as characteristic 
examples chains with one steel and two aluminum beads, as well as ones with 1 steel and three alu- 
minum beads, although we briefly touch upon other configurations as well, illustrating their similar 
characteristics. We analyze the higher order gaps that emerge in such settings and explore the 
intrinsic localized modes that bifurcate from the edge of the upper bands. A generic surprising fea- 
ture of such states is that they appear to be more robust than their counterparts bifurcating from 
the edges of the lower bands. Direct numerical simulations, using driving of the system at suitable 
frequencies through an actuator or taking advantage of the modulational instabilities of extended 
band edge states in the system illustrate the spontaneous formation of localized modes within the 
corresponding nearest gaps. In these numerical experiments, we also account for the presence of 
dissipation and analyze its potential role. 



INTRODUCTION 

The theme of "granular crystals," i.e., of chains 
composed of elastic (typically spherical) particles in 
close contact with one another, is one that has been 
gaining considerable momentum over the past few 
years [l], Q ■ Of particular interest in this setting is 
the nonlinear interaction between the grains, which 
features a so-called Hertzian contact between the ad- 
jacent beads in the chain upon compression and a 
zero tensile strength. The Hertzian interaction po- 
tential depends on the relative displacement of the 
bead centers raised to the 5/2 power for spherical 
beads. One of the remarkable features of this sys- 
tem is that it can be tuned to possess a dynamic 
response ranging (depending on the static compres- 
sion applied to the chain) from weakly all the way 
to strongly nonlinear [l|, Q . It is precisely this tun- 
ability and straightforward experimental implemen- 
tation of such chains that enables a wide range of 
nonlinear dynamical studies in them, featuring the 
interplay of discreteness, nonlinearity and period- 
icity (or lack thereof). Additionally, such systems 
have shown considerable potential for a wide range of 
engineering applications including shock and energy 
absorbing layers , actuating devices Q , acoustic 
lens es Hi , acoustic diodes Q , and sound scramblers 
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Our emphasis in the present work will be on the 
coherent nonlinear waveforms that can trap the en- 
ergy in such systems. Among the most prominent 
ones such that have been reported in the literature 
are the solitary waves 0, 0, 12, EH, shock waves 
14 , ]j| and very recently discrete breathers 1||- 18 1 



(see also related discussions in Refs. [19|, |20(). Our 
emphasis in the present work will be on the discrete 
breathers, which have a time honored history in a di- 
verse host of settings. Such intrinsic localized modes 
constitute generic excitations of nonlinear dynamical 
lattices, which are (typically) exponentially localized 
in space, while they oscillate periodically in time. 
The discrete breathers appear to be relevant to a 
wide array of thematic areas spanning superconduct- 
ing Josephson junctions 



2l| . photonic crystals 



biopolymersj2 
romagnets [2 
rays [26[, among others [27 



, charge-transfer solids [24j , antifer- 
and micromechanical cantilever ar- 



It has been argued 2, 17 1 that discrete breathers 
arc likely to be absent in monoatomic granular 
chains 281. Hence, their first theoretical investi- 



gations and experimental implementations in the 
above works took place in the context of heteroge- 
neous systems in the form of dimers. Such arrays 
consisting of glued (29[, welde d l30l . and elastically 
compressed spherical particles 1^ 31-3^1 have been 
shown to exhibit tunable frequency vibrational band 
gaps. The works of Refs. [16|, |l7| took advantage 
of the modulational instability [27[ of the mode at 
the bottom edge of the optical band in this type 
of chain to produce dynamical localization in the 
form of discrete breathers whose frequency of oscilla- 
tion resided in the nearest bandgap (while the length 
scale of their exponential localization depended on 
how deep in the bandgap their frequency was lying) . 

However, in recent experimental works [34T - [36j ]. 
experimental accessibility has been demonstrated to- 
wards states with generalized periodicities that are 
not restricted merely to dimers, but may include 



trimers (such as so-called 2 : 1 configurations, with 
2 beads of one material and one bead of another, 
or 1 : 1 : 1 chains consisting of 3 different material 
types) or even quadrimers (with 3:1, 1:1:1:1, 
2 : 2, 2 : 1 : 1 constituting possible combinations 
of materials therein). While solitary traveling waves 
have been identified in such chains in the work of 
Refs. (Hill, and their linear spectrum under pre- 
compression has been illustrated to contain multiple 
bandgaps in the recent work of Ref. (36j . nonlin- 
ear breathing excitations have not, to the best of 
our knowledge, been identified therein. This consti- 
tutes the focus of the present study, which is partly 
motivated from the analog of such investigations in 
nonlinear Schrodinger models with optical lattices 
(related to Bose-Einstein condensates) whereby gap 
solitons in each of the gaps of the optical lattice po- 
tential have been identified in the recent work of Ref. 
371. 



Herein, upon presenting the general theoretical 
setup (Sec. 2), we restrict ourselves to arguably 
the simplest trimer and quadrimer configurations, 
namely, the experimentally accessible 2 : 1 chains 
with two aluminum and one steel bead (Sec. 3), 
as well as their 3 : 1 counterparts with three alu- 
minum and one steel bead (in Sec. 4). Nevertheless, 
the generality of our conclusions is shown to be sup- 
ported by the examination of additional configura- 
tions, such as two steel beads and one aluminum, or 
a trimer that involves aluminum, steel, and brass. 

We commence their examination by analyzing the 
corresponding linear spectrum. We find here (as well 
as more generally) that configurations with a peri- 
odicity of P beads in such granular systems possess 
P — 1 finite gaps and a semi-infinite one (extend- 
ing from the top of the uppermost optical band to 
oo). Within the resulting 2 finite gaps of the trimer 
chain and the 3 of the quadrimer one we seek to 
identify discrete breather states and systematically 
continue them within the relevant bandgap. We 
believe that this is the first demonstration of such 
higher bandgap states in granular systems. More- 
over, such states have a rather surprising character- 
istic: They appear to be considerably more robust 
than their lower gap counterparts; that is, they have 
wider parametric intervals of stability and even when 
unstable their instability growth rates are typically 
considerably lower than their lower bandgap coun- 
terparts. Finally, in each case, to demonstrate the 
experimental realizability of the pertinent states, we 
report dynamical simulations under experimentally 
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FIG. 1: (Color online) A schematic of our 2 : 1 
aluminum-steel configuration. 



realistic conditions; in that light we also account for 
the presence of dissipation and elucidate its poten- 
tial role. Wc induce the formation of these nonlinear 
gap solitary wave (or gap breather) states through 
the excitation of a modulationally unstable mode of 
the system at the frequency of the bottom of the 
optical band edge (different periodicities exhibit dif- 
ferent numbers of such band edge modes). This is 
either done by initializing the system at the spatial 
profile of the extended mode of the proper frequency, 
or, more realistically, by driving ("actuating") the 
system at the relevant frequency; cf. also the exper- 
imental realization of the lowest bandgap breathers 
in Ref. flij ]. Finally, in Sec. 5, wc summarize our 
findings and present our conclusions. 



THEORETICAL SETUP 

We will be focusing on a one-dimensional (ID) 
granular crystal containing steel and aluminum 
spheres compressed uniaxially by a force Fq = 
22. 4N (this ensures the existence of a linear spec- 
trum, which is critical for enabling the generation of 
breathers within the spectral gaps). This granular 
crystal is modeled by a system of coupled nonlinear 
oscillators; the equation of motion for the displace- 
ment u n of the nth bead is given by [l], Q 

m n u n = A n -i tn [S n - ltn - (u n - ii„-i)]f 2 

- A n .n+l[Sn,n+l ~ ("n+1 ~ Ufi)]+ 2 , (1) 

where m„ is the mass of the nth bead, A n ^ n+ i is the 
coefficient of elasticity between beads n and n + 1, 
and <J n ,n+i — (-PoA4-ri,n+i) 2 ^ 3 is the amount of static 
overlap of the beads n and n + 1 (due to the precom- 
pression) under equilibrium conditions. The coeffi- 
cient A niTl+ i depends on the material and geometric 
properties of the beads n and n + 1; for sphere-to- 
sphere contact, which will be the case of interest 
herein, the coefficient A njn+ i is given by 



4E n E n+ i 



r n r n +i 
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An > n+1 3E n+1 (l vl) + 3E n (l 



(2) 
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where r n , E n , and v n denote, respectively, the ra- 
dius, elastic (Young's) modulus, and Poisson ratio 
of the nth bead. The bracket [s]+ of Eq. (Q} is the 
positive part of s, that is, 

[s] + = max(0, s) 

This accounts for the fact that if two beads are not in 
contact, they do not exert any force on each other 
(zero-tension). In addition, each bead (away from 
the boundary) has a local energy density associated 
with it, given by 

e n = ^rn n Un + -[V n - 1 (u n -u n _i) + V n (u n+1 -u n )], 

(3) 

where the interaction potential V n (u n +i — u n ) is 
given by 

V n (u n+ 1 - U n ) =-j4„ n+1 [A nn+ i - (u n +l - Un)]+ 2 
5 

- jU„,„ +1 A^ +1 . (4) 

We note here that the sum of ([3]) over all the beads 
gives us the Hamiltonian of the system, and the 
equations of motion can then be derived in the usual 
way (using Hamilton's equations). In Table HI we 
present an overview of the material properties used, 
in sync with earlier experimental results [lij ]. The 
trimer lattice will have repeating unit cells contain- 
ing a 'aluminum-aluminum-steel' (a-a-s) pattern; see 
Fig. [T] We will also consider chains with four 
beads per unit cell, i.e., an 'aluminum-aluminum- 
aluminum-steer (a-a-a-s) pattern. Although addi- 
tional heterogeneous configurations are possible, per 
our discussion in the previous section, our results in- 
dicate that these settings accurately reflect the po- 
tential for classes of additional gap breather modes. 

A-A-S CHAIN 
Linear Spectrum 

Except for some degenerate (non-generic) cases, 
which depend on the materials used, the periodic- 
ity of the lattice determines the number of distinct 
bands in the dispersion relation; each band is a fam- 
ily of solutions to a polynomial equation (specifically 
a P x P matrix eigenvalue problem) , and the degree 
of the polynomial is the same as the period of the 



Material Steel Aluminum 



Radius (mm) 


9.53 


9.53 


Density (kg/m 3 ) 


8027.17 


2700 


Elastic modulus (Pa) 


2 x 10 11 


7 x 10 10 


Poisson ratio 


0.3 


0.35 



TABLE I: The material properties used in our 
simulations; cf. the experimental parameters used 
inRef. 
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FIG. 2: (Color online) Left: The linear spectrum 

for a chain of 80 beads with period 3 
(aluminum-aluminum-steel). There are three linear 
bands separated by two finite gaps. Right: The 

eigenmodes from the bottom edges of {top to 
bottom) the second and first optical bands of the 
chain. 



lattice. Thus, in general, the number of finite gaps 
in the spectrum is one less than the period of the 
lattice (there is generically also a semi-infinite gap 
extending from the top of the highest optical band to 
+00). Figure |2] shows the spectrum for an a-a-s lat- 
tice along with sample eigenmodes from the bottom 
edges of the bands; these will be used as the initial 
seeds for our nonlinear solvers that will be used to 
identify the discrete breather modes in the gaps of 
the corresponding linear spectra. 

Some of the features worth mentioning within the 
relevant spectrum are the following. Firstly, we 
can observe that the free boundary conditions used, 
along with lighter masses at each end of the chain, 
induce the presence of some localized modes centered 
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around the edge of the chain with frequencies within 
the linear gap. Such "surface modes" have been dis- 
cussed earlier in the context of Refs. 17, 18| and will 
not be further elaborated herein. The modes at the 
bottom edges of both optical bands are reminiscent 
of the ones reported in the dimer case in Ref. [l6| 
in that the "heavy" steel beads are immobile while 
the "lighter" aluminum ones are the ones oscillating. 
As may be anticipated, in the lower band, these alu- 
minum pairs oscillate in phase, while in the upper 
one, they oscillate out of phase. It should be noted 
that a detailed analysis (through explicit dispersion 
relations) of the linear spectrum in the case of an 
a-a-s lattice can be found in the recent theoretical 
and experimental work of Ref. [36| ; hence we refer 
the interested reader there for further details on the 
linear picture. 



Existence and Stability of Discrete Breathers 

Using an initial guess from the bottom of the band 
edge, we performed a numerical single-parameter 
continuation of a Poincare map through each spec- 
tral gap using frequency as our continuation param- 
eter. See the Appendix for more details about this 
method of continuation, and the Floquet multipli- 
ers. Sample breather profiles from these continua- 
tions are shown in Fig. [3j It can be observed that 
similarly to the profiles of the linear modes from 
which they stem, the breather modes bear an ap- 
parent in-phase excitation of their aluminum beads 
in the lower gap, and an out-of-phase one in the up- 
per gap. In these nonlinear modes, the heavier steel 
beads gradually become excited as well. Further- 
more, the deeper into the respective gaps the modes 
get (in terms of frequency) , the more localized they 
become, as is evident by the strongly localized form 
of the gap breather profiles for the chosen frequen- 
cies lying around the middle of the respective gaps. 

In parallel to the numerical continuation, we 
also examined the linear stability of each obtained 
breather solution via Floquet analysis. The rele- 
vant results are shown in the full linearization spec- 
tra (corresponding to the breather profiles) of Fig. 
[3] and the continuation results of Fig. |4j When- 
ever any of these Floquet multipliers A^ satisfies 
| Ai| > 1, the corresponding breather is linearly un- 
stable. There are two types of instabilities that we 
observe in this Hamiltonian system. The one that 
stems from a purely real Floquet multiplier with 



magnitude greater than one (which forms a pair with 
another such with magnitude less than one) , we will 
call a real instability. On the other hand, there also 
exists the possibility of a quartet of Floquet multi- 
pliers outside the unit circle. These modes bearing a 
nonzero imaginary part are typically weaker, and we 
will call them oscillatory instabilities, as they are as- 
sociated with complex eigenvalues and exhibit both 
oscillatory and growth dynamics. 

As has been observed previously in Refs. []~6l — 
HH, a change in the monotonicity of the energy as 
a function of frequency is directly associated with a 
real instability for the case of gap breathers. Such 
an instability was earlier observed for the first gap 
in the dimer case in Refs. [H, |17| . and a similar 
finding arises in the lower finite gap of the trimer 
case considered herein. Furthermore, the structure 
of the oscillatory instabilities appears to be similar 
between the two settings. Interestingly, once the fre- 
quency becomes lower than a critical point, call it f c , 
sharp resonances with extended modes emerge sim- 
ilar to what was observed in Ref. [13] and the gap 
breathers become highly unstable (hence, we do not 
consider such modes for lower frequencies herein). 
This critical point / c , as it turns out, is the point at 
which the 2nd harmonic of the breather frequency 
lies in the highest linear band, and thus the breather 
is resonating with the linear modes at the second 
harmonic. This is confirmed by observing that the 
breathers found with frequency lower than f c have 
excited tails, which oscillate with a frequency twice 
that of the breather. 

However, the key observation of the present work 
is that the modes of the higher gap have a signifi- 
cantly narrower frequency range for which they ex- 
hibit the real instability (in fact, for different mate- 
rial properties the instability may be entirely absent; 
see below). Furthermore, for the oscillatory insta- 
bilities too, the growth rates appear to be consid- 
erably weaker than their lower gap analogs in Fig. 
[31 apparently leading to considerably more robust 
higher gap states. This is a rather surprising obser- 
vation and may perhaps be attributed to the more 
pronounced out-of-phase structure of the higher gap 
modes, which is well known in other classes of dy- 
namical lattice models to offer enhanced stability as 
well [H]. 

For completeness, we also performed a continua- 
tion through the spectral gap, starting with an initial 
seed of the surface mode in the gap, which is present 
because of the free boundary conditions, and the 
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FIG. 3: (Color online) Top: A sample breather 
profile from the upper band gap, along with its 
Floquet multipliers. Bottom: A sample breather 
from the lower band gap, along with its Floquet 
multipliers. 

light masses on the ends. These surface breathers 
were incredibly robust along the entire continua- 
tion, with all Floquet multipliers staying within a 
0.01 radius of the unit circle. The energy of these 
breathers increased until it closes in on the branch 
from the bulk of the chain, and then follows closely 
that branch through the rest of the continuation. 
The profile of the surface breather corresponds al- 
most exactly to the profile of the corresponding bulk 
breather, as shown in Fig. This suggests that sur- 
face breathers are effectively bulk breathers which 
have been shifted to the boundary. 
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FIG. 4: (Color online) Top: The magnitude of each 
Floquet multiplier is shown as a function of its 
frequency. Purely real Floquet multipliers are 
shown in black, while Floquet multipliers 
associated with nonzero imaginary part are shown 
in gray. Bottom: The energy of each solution as a 

function of frequency. The sharp corner in the 
energy happens at the critical frequency f c , where 

the breathers start resonating with the linear 
modes due to the second harmonic. The horizontal 
scale is the same for both panels, and the frequency 
ranges of the linear bands are shaded for clarity. 
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FIG. 5: (Color online) The surface breather with 
frequency of 7.9764 kHz is shown shifted and 
overlaid on the bulk breather of the same frequency. 



Dynamics and Modulational Instability 

We now turn to the examination of the potential 
experimental feasibility of exciting such nonlinear lo- 
calized modes in the 2:1 chain. To that effect, we will 
examine a number of direct numerical simulations of 
the chain evolution under suitable initial conditions 
or external drive. In what follows, we focus, in par- 
ticular, on the breathers from the upper band gap, 
since they constitute the most robust (and, hence, 
most likely to be experimentally observable) struc- 
tures of the system. 

First, using the (properly rescaled) eigenmode 
from the bottom edge of the highest band as an ini- 
tial condition, we observed deviations from the opti- 
cal eigenmode which were reinforced by the nonlin- 
earity, leading to the eventual breakup of the wave- 
form. As a result of this process, called modula- 



tional instability, a breather localized near the cen- 
ter of the lattice was formed. This is shown in the 
top panel of Fig. [6j which illustrates just how ro- 
bust these breathers are. A power spectral density 
(PSD) analysis of these breathers revealed that they 
did indeed acquire frequencies in the upper band- 
gap; in this particular case, the dynamically result- 
ing breather had a frequency of about 8.23 kHz (4l| . 
The main conclusion from this and other similar 
evolution runs that we performed is that the top- 
band-edge breathers are prone to the manifestation 
of the modulational instability in ways similar to 
their lower gap analo g in the dimer case experiments 
and analysis of Ref . [16( ■ 

Since dissipation is ubiquitous in experimental se- 
tups (39l | . we then added linear on-site damping to 
each bead in order to determine whether dissipa- 
tion would have any effect on the formation of these 
breathers due to modulational instability. With this 
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FIG. 6: (Color online) Top: Modulational 
instability of the first eigenmode of the upper most 

band, in a chain with 'aluminum - aluminum - 
steel' cells. A robust, long-lived breather is formed. 
Bottom: The same as the top panel, except in this 
case, dissipation is turned on, with r = 2 ms. No 
breather forms. 

dissipation term, the equations of motion arc modi- 
fied to be 



3/2 

+ 

A„.„._|_1 I5n.-n.4--l — {u n -\-l ~ u n)]-l 



m n u n = A n _i jn [8 n -i, n - (u n - u„-i)] 

"tn+l 



in,n+l [^n,n-t-l 



(5) 



where r = 2 ms was chosen to best simulate the ef- 
fects of dissipation in an experimental setting jioj . 
The experiments above involving the linear mode 
as an initial condition were redone with dissipa- 
tion turned on, and the results were strikingly dif- 
ferent. Instead of observing modulational instabil- 
ity and the formation of a breather, the oscillations 
simply decayed exponentially. It appears that the 
timescales over which the solutions decay are too 
short to observe modulational instability. This com- 
parison of the dissipative and Hamiltonian case can 
be seen between the bottom and top panels of Fig. 
|6j It is worth noting that in other simulations with 
much weaker dissipation, we did observe the forma- 
tion of breathers, which were, in fact, more robust 
than those of the purely Hamiltonian case. This 
may become especially relevant if experimental tech- 
niques are devised that decrease the effects of dissi- 
pation. See also below for a detailed discussion on 
the effects of dissipation for an experimentally rele- 
vant technique of producing such discrete breathers. 

Lastly, we investigated the possibility of nucleat- 
ing one of these higher gap breathers using an actu- 
ator, a technique that has been experimentally em- 



ployed in the dimer case in Ref. [lfj. Rendering 
the model more realistic, by keeping the dissipative 
term, we added a periodic (in time) drive beside the 
first bead in the chain (recall that previously, the 
chain had free boundary conditions). We set this 
"actuator," a harmonically oscillating bead at the 
left (n = 0) boundary, to ramp up linearly to its full 
amplitude, then remain on for a period of time, and 
finally shut off, reverting to a fixed boundary. 



a [Jr) sin(uj d t) for < t < T r 
\asin(uj d t) for T r < t < T d 

for t > T d 



(6) 



where a is the actuator amplitude, T r is the ramp- 
up time, c — 1 for linear ramping, uj d is the actuator 
frequency, and T d is the actuator duration. The ac- 
companying intuition here is that with the actuator 
pumping at a frequency right at the bottom edge of 
the upper band, it is possible to excite, via modula- 
tional instability, the formation of a breather some- 
where in the chain. With extensive experimenta- 
tion, we found that the only breathers formed were 
on the surface of the chain, but as we mentioned 
above, these surface breathers do have similar pro- 
files to those of the bulk breathers. Figure [7] shows 
the results of one such experiment. Here, the actua- 
tor ramp time was 20 ms, the maximum amplitude 
of the actuator was 9 x 10 -7 m, and the actuator 
was allowed to run for 80 ms before being shut off. 
A breather was formed right near the surface, which 
extended about two unit cells into the chain. This 
breather had a frequency of about 8.20 kHz, which 
is well into the upper band gap. The profile of this 
breather is the same, at least qualitatively, with the 
corresponding exact solution found by the continu- 
ation at this frequency. 

Comparison to Other Trimers 

In this subsection the results of the a-a-s chain 
are compared to the results for two other trimers: 
one consisting of an 'aluminum-steel-steel' (a-s-s) 
pattern, and the other consisting of an 'aluminum- 
steel-brass' (a-s-b) pattern. In the results of each 
of these continuations, shown in Figs. [8] and [9l we 
see that for each of these configurations, the proper- 
ties of the breathers are quite similar. More specif- 
ically, we identify as before three pass bands and 
between them two finite gaps (as well as a semi- 
infinite gap above the highest band) in the linear 
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FIG. 7: (Color online) Top: The velocities of an 
experiment with an actuator having frequency right 

at the lower edge of the third band (fd — 9.084 
kHz), amplitude 9 x 10 m, and dissipation (with 
r = 2 ms). Middle: The profile of the chain at the 
time indicated by the vertical bar in the upper 

panel. This profile corresponds, at least 
qualitatively with both the surface breather and 
bulk breather of the same frequency. Bottom: The 
PSD of the dynamical force experienced by bead 3. 
The linear bands are colored in so that the gaps are 

visible. There is clear evidence of a frequency 
within the linear band gap, and the breather has a 
frequency of approximately 8.20 kHz. 



spectrum. It is within these finite gaps that we fo- 
cus on, regarding the potential existence of discrete 
breathers. The lower gap has breathers that are 
mainly unstable, and not very localized, while the 
upper gap has breathers that are only weakly unsta- 
ble, and are highly localized. The main difference 
between these two configurations and the previous 
(s-a-a) configuration is the absence of any real in- 
stability in the breathers of the higher gap, through 
a pair of Floquet multipliers exiting from unity onto 
the real axis,. i.e., the breather energy is found to 
have a monotonic dependence on the corresponding 
frequency. Additionally, that critical frequency f c 
at which resonances with the second harmonic ap- 
pear, due to the differences in the linear spectrum, 
is shifted in both cases. In particular, for the a-s- 
b chain, f c is shifted almost to the bottom edge of 
the second band, which prevents the numerical con- 
tinuation from finding more than a small segment 
of breathers. Nevertheless, we highlight that all the 
principal characteristics of both the linear spectrum 
and of its gaps being conducive to the formation of 
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FIG. 8: (Color online) The results for the 
continuation of breathers in an 
'aluminum-steel-steel' chain. Top left: The energy 
of each breather solution is plotted as a function of 
frequency. Bottom left: The complex modulus of 
each Floquet multiplier is shown as a function of 
its frequency. Top right: The linear spectrum of 
the a-s-s chain is shown. Bottom right: Two 
sample profiles of breathers in this chain are shown, 
one from each gap. 



the nonlinear waveforms of interest (with the same 
general structural features as before) are shared be- 
tween these alternative configurations and the orig- 
inally examined a-a-s one, attesting to the generic 
nature of our conclusions. 



The Role of Dissipation 

In this subsection, we examine the role of dissipa- 
tion on the results associated with the actuation of 
the chain at its boundaries. In particular, we per- 
form a series of experiments which are identical to 
those of Fig. [71 but now the dissipation is made ei- 
ther 10 times larger, i.e., t = 0.2ms, as shown in 
Fig. [El or it is 10 times smaller, i.e., r = 20ms, as 
presented in Fig. [IT] In the former case, it is clear 
that the dissipation completely overwhelms any po- 
tential nonlinear structure formation, coming into 
play at a very early time scale. It is for that reason 
that the dynamics of the fifth site is purely oscilla- 
tory, and the relevant oscillation frequency, as illus- 
trated in the corresponding power spectral density, 
is found to be the driving frequency at the bottom 
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FIG. 9: (Color online) The results for the 
continuation of breathers in an 
'aluminum-steel-brass' chain. Top left: the energy 
of each breather solution is plotted as a function of 
frequency. Bottom left: the complex modulus of 
each Floquct multiplier is shown as a function of 
its frequency. Top right: the linear spectrum of the 
a-s-b chain is shown. Bottom right: Two sample 
profiles of breathers in this chain are shown, one 
from each gap. 



edge of the third band. On the contrary, and in a 
way more reminiscent of the dynamics of Fig. [Jj the 
dynamics of the case with r = 20 ms is more con- 
ducive to the emergence of additional cigenfrcquen- 
cies through the modulational instability mechanism 
highlighted previously. It is for that reason that the 
evolution of the force at the fifth site appears to in- 
volve multiple frequencies, as is confirmed also by 
the corresponding power spectral density. Further- 
more, as expected, in addition to the driving fre- 
quency, there exists a maximally excited frequency 
within the second gap which is associated with the 
discrete surface breather emerging in the dynamics 
of the top panel. 

In conclusion, for the range of dissipations that 
have been recently used to match experimental re- 
sults [40( | . as well as for weaker dissipations, the 
waveforms identified herein are expected to persist. 
On the other hand, for the case of materials with sig- 
nificantly more dissipative dynamics (and possibly 
effects of plasticity that additionally drain the en- 
ergy pumped into the system) , these coherent struc- 
tures cannot develop. 
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FIG. 10: (Color online) The diagnostics shown here 
are the same as in Fig. [TJbut now for a dissipation 
which is 10 times larger, i.e., r = 0.2 ms. 
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FIG. 11: (Color online) The diagnostics shown here 
are the same as in Fig. [TJbut now for a dissipation 
which is 10 times smaller, i.e., r = 20 ms in this 
case. 



A-A-A-S CHAIN 
Linear Spectrum 

In this section we present similar results to the 
above a-a-s chain, but for a chain with periodicity 
involving P = 4 sites, with each unit cell following 
an a-a-a-s pattern. Similarly to the a-a-s case, we 
start by illustrating in Fig. [T^] the linear spectrum 
of a chain with 160 beads, which has four bands and 
three finite gaps. Examples of the linear modes at 
the bottom of each of the bands (second, third, and 
fourth, respectively) are shown in the right panel of 
Fig. HU It can be seen that at the bottom of the 
second band, the mode involves an in-phase excita- 
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FIG. 12: (Color online) The linear spectrum for a 
chain with period 4 
(aluminum-aluminum-aluminum-steel) . The 
spectrum features three finite gaps separating the 
four existing pass bands. Within the gaps, the 
(free) boundary conditions induce a number of 
localized surface modes (at the chain boundaries). 
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FIG. 13: (Color online) Top: A sample breather 
profile from the third band gap, along with its 
Floquet multipliers. Middle: A sample breather 
profile from the middle (finite) band gap. Bottom: 
A sample breather profile from the lowest band 
gap. 



tion of the aluminum beads (while the steel ones are 
nearly stationary). At the bottom of the third band 
every other bead (one steel and one aluminum) are 
essentially stationary, while the remaining two alu- 
minum beads that straddle the stationary one are 
out of phase with respect to each other. Finally 
at the bottom of fourth band, each aluminum bead 
is essentially out of phase with its neighboring one 
(while the heavier steel beads remain stationary). 



Existence and Stability of Discrete Breathers 

Once again, a continuation through each gap 
was performed, starting from the corresponding 
band edge modes and, accordingly, three families of 
breathers were discovered. A sample breather from 
each gap is shown in Fig. [131 It can be seen that 
each of the breathers preserves a structure reminis- 
cent of the corresponding linear mode from which it 
stemmed, while forming an increasing displacement 
step due to the asymmetry of the potential (as dis- 



cussed earlier in Ref. [H(). As before, the breathers 



acquire an increasing degree of localization as the 
frequency gets deeper into the corresponding band 
gap. 

Finally, the linear stability of each breather was 
computed via Floquet analysis, and the magnitude 
of the Floquet multipliers was plotted as a func- 
tion of frequency in Fig. 1141 Remarkably, once 
again breathers in the uppermost finite gap are 
the most robust among the three families. In par- 



ticular, they feature no real instability whatsoever 
and, in fact, have an interval of complete stabil- 
ity, while their oscillatory instabilities are only very 
weak when present. On the other hand, the middle 
gap breathers are highly unstable, featuring a strong 
real instability in a large fraction of the correspond- 
ing band gap (and even a second real instability close 
to the top of the second band). Nevertheless, they 
also possess an interval of stability in the immediate 
vicinity of of the bottom of the third band; there 
is also a mild real instability in that neighborhood 
associated with an energy slope change; see Fig. fT4l 
and cf. Ref. [l7j]. Finally, within the lower gap, the 
discrete granular breathers are so spatially extended 
that boundaries seem to interfere with our numeri- 
cal continuation demanding a very large domain to 
accurately perform the relevant computation. Nev- 
ertheless, these states also appear to be subject to 
strong real instabilities. Thus, once again the top 
finite gap breathers appear to be the most robust 
nonlinear entities of the system. 



Dynamics and Modulational Instability 

We also investigated the potential for exciting dis- 
crete breathers in the 3:1 chain. Once again, we used 
a rescaled linear mode from the bottom edge of each 
optical band as an initial condition in our numerical 
integrator, and observed that modulational instabil- 
ity did indeed form a breather localized near the 
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FIG. 14: (Color online) Top: The magnitude of 
each Floquet multiplier is shown as a function of 
its frequency. Purely real Floquet multipliers are 
shown in black, while Floquet multipliers with 
nonzero imaginary part are shown in gray. Bottom: 
The energy of each solution as a function of 
frequency. The horizontal scale is the same for 
both panels. In both plots, the frequency ranges of 
the linear bands are shaded for clarity. 
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FIG. 15: (Color online) The results of a numerical 

experiment of the modulational instability of a 
linear mode at the bottom edge of the third (top) 
and second (bottom) optical band. In each case a 
breather was formed, which was localized near the 

center of the lattice. These breathers had 
frequencies of approximately 9.12 kHz and 7 kHz, 
respectively, in the second and third finite band 
gaps. 



center of the chain. Figure [TS] shows the results of 
these experiments. A PSD analysis shows that the 
frequency of each breather is located well into the 
corresponding band gaps. 

In order to excite the relevant modes, we have once 
again resorted to the actuation of the band edge of 
(the bottom of) the different pass bands, in line with 
our description in the previous section (we similarly 
use a ramp and eventually cut off the driving of the 
chain boundary). We have also included dissipation 
in the model to render our results closer to the exper- 
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FIG. 16: (Color online) Top: The velocities of a 
numerical experiment with an actuator having 
frequency at the bottom of the fourth band (9.45 
kHz), amplitude 5 x 1.5~ 6 , and dissipation (with 
t = 2). The chain in this case had 159 beads in 
order for boundary reflection effects to be lessened, 
and a clearer picture obtained. Bottom: The PSD 
of the force time series of bead 5 is shown. The 
linear bands are colored in so that the gaps are 
visible. There is clear evidence of a frequency 
within the top band gap, and the breather has a 
frequency of approximately 8.89 kHz. 



imentally anticipated ones. Our findings arc shown 
in Figs. H6l and[T7 ] respectively, for the top and mid- 
dle finite gaps. The conclusion in both cases is, in 
fact, similar, namely the drive of the mode at the 
bottom of the optical band edge induces a modula- 
tionally unstable evolution, in turn, leading to local- 
ization within the respective band gap. As a result 
the former figure features a breather of frequency 
8.89 kHz, as a result of the dynamical evolution, 
while the latter one similarly possesses a localized 
mode of frequency 7.03KHz within the middle gap. 
These results indicate that the experimental realiza- 
tion of such higher gap breathers should certainly be 
within reach with currently accessible experimental 
settings. 



CONCLUSIONS AND FUTURE 
CHALLENGES 

In the present work, we have extended earlier con- 
siderations focusing on precompressed elastic dimers 
(featuring nonlinear discrete breather modes within 
their single finite spectral gap). In particular, we 
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FIG. 17: (Color online) Top: The velocities of an 
experiment with an actuator having frequency at 
the bottom of the third band (7.59 kHz), amplitude 

8 x 10 -7 , and dissipation (with t = 2). Bottom: 
The PSD of the force time series of bead 5 is shown. 
The linear bands are colored in so that the gaps are 

visible. There is clear evidence of a frequency 
within the middle finite band gap, and the breather 
has a frequency of approximately 7.03 kHz. 



hand, from the theoretical perspective, it would be of 
particular interest to generalize the present consid- 
erations not only to alternative one-dimensional het- 
erogeneous configurations, but perhaps more appcal- 
ingly to the substantially more challenging setting of 
two-dimensional granular crystals. Understanding 
the bandgaps of such two-dimensional chains and 
the associated intrinsically localized modes thereof 
would be a theme of particular relevance for theo- 
retical and experimental studies alike. 
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granular chain setup and its modeling. 



Appendix: Poincare Continuation 



have considered the linear and associated nonlinear 
properties of, arguably, one of the simplest forms 
of trimer and of quadrimcr chains, namely a 2:1 
aluminum-steel chain and a 3:1 such chain. We 
have illustrated the presence of 2 finite band gaps 
in the first case and of three such in the second. 
Within these bandgaps, we have been able to iden- 
tify (starting from the bottom edge of the upper 
band in each case) a sequence of gap breather local- 
ized modes with increased localization properties, as 
the frequency of their temporal oscillation proceeded 
deeper into the linear gap. We have also showcased 
the highly nontrivial feature that the states of the 
higher gaps appear, in fact, to be more robust than 
the ones of the lower gaps, presumably a feature 
that should be attributed to their out-of-phase spa- 
tial structure. Finally, we have also demonstrated 
the spontaneous emergence of such states, as a re- 
sult not only of the modulational instability of the 
extended modes of the bottom edge of the optical 
bands of the linear spectrum, but also as a result of 
the experimentally realizable scenario of the system 
being driven at such band edge frequencies. 

We believe that these findings pave the way to- 
wards the experimental realization and observation 
of such states and more generally towards a system- 
atic classification of gap breather modes in elastic 
systems with multiple band gaps. On the other 



In order to find periodic solutions to this system, 
we are searching for solutions u so that u n (0) = 
UniT) where T is the fixed period of the solu- 
tion. This gives us a Poincare map, P(uq) = 
u(0; uo) — u(T; uo), where uq is the initial condition, 
and u(t; uq) is a solution to system of ODEs at time 
t with initial condition uq. A periodic solution with 
period T (frequency 1/T) will be a root of P. Since 
any point on a periodic solution will also be a root, 
we add the additional constraints that u n = in 
order to "pin" the solution to a point in the do- 
main of the Poincare map. To calculate the numeri- 
cally exact (up to a set relative error, usually 10 -10 ) 
profile of each periodic solution, we used a Newton- 
Raphson solver to calculate the fixed-points of P. 
The Jacobian for the Newton's method is given by 
/ — V(T), where V(T) is the monodromy matrix, 
which we obtain by integrating the variational equa- 
tion 

j t V = J(u,u)V, 

simultaneously with the equations of motion ((T|), 
where J is the Jacobian of these nonlinear equations, 
and 1^(0) = / (i.e., the intial condition of the mon- 
odromy matrix is the identity matrix). The eigen- 
values of the monodromy matrix also give us the 
Floquet multipliers of the solution, which convey in- 
formation about the linear stability of the obtained 
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periodic orbit. Floquet multipliers outside the unit 
circle give rise to instabilities, while Floquet multi- 
pliers on the unit circle lead to marginal stability, 
and Floquet multipliers purely inside the unit circle 
give rise to asymptotic stability. 

Once we have our numerically exact solution u, 
which is hopefully a localized DB, we then change 
the frequency (and thus the period) of our desired 
solution, and using u as an initial guess, we again 
use the Newton solver to find another exact solution. 
An entire family of solutions can then be found by 
'continuing' this solution through the frequency gap. 
For a review of the many methods of finding DBs, 



see Ref. 271 and the references therein. For more 



details about numerical continuation, see Ref. 42 
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